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In this report, progress in the investigation of empennage 
buffeting is reviewed. In summary, the following tasks have been 
accomplished : 

(1) Relevant literatures have been reviewed. 

(2) Equations for calculating structural response have been 
formulated. 

(3) Root-mean-square values of root bending moment for a 65- 
degree rigid delta wing have been calculated and compared 
with data. 

(4) Water- tunnel test program for an F-18 model has been 
completed. 

Items (1) - (3) are described in more detail in Appendix A, while 
item (4) is presented in Appendix B. 




Appendix A: 

Investigation of Buffeting by LEX Vortices 


by 


C. Edward Lan and I. G. Lee 
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INTRODUCTION 


Buffeting flow arises when flow separation occurs on an 
airplane. The resulting flow field is highly turbulent, thus 
producing fluctuating pressures on lifting surfaces in the detached 
flow region. Boundary layer separation is perhaps the most common 
source producing buffet on most conventional configurations. 

Research in this area has been quite extensive and involved 
measurements of fluctuating pressures on models together with some 
theoretical methods to extrapolate these results to full-scale 
vehicles (see, for example. Refs. 1-3). Frequently, these pressure 
measurements are made on a conventional "rigid" model, instead of an 
aeroelastic one, because the latter can not withstand high enough 
dynamic pressures to be realistic. Based on this consideration, 
several theoretical methods to use these pressure measurements to 
predict buffet response have been developed. Some of these methods 
will be reviewed later. Review of some test results can be found in 
References 4 and 5? and of theoretical methods, in References 6 and 
7. 

Of particular interest in the present investigation is the 
buffeting caused by leading-edge vortices on slender wings. Test 
results showed that 

(a) buffeting was low before vortex breakdown and became 
severe after that (Ref. 8 and 9); 

(b) high-frequency buffeting was caused by boundary layer 
fluctuation, and leading-edge vortices produced mainly 
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low-frequency fluctuation (Ref. 10); 

(c) the results were not sensitive to Reynolds numbers (Refs. 
10 and 11), so that flight and tunnel measurements could 
be well correlated (Ref. 12); 

(d) buffeting at vortex breakdown was associated with the wing 
response at the fundamental mode (Ref. 8). 

One conclusion . from this early-day research on leading-edge 
vortces was that the buffeting induced by vortex breakdown would 
mostly be academic because a slender-wing airplane would normally 
not operate in the vortex-breakdown region of angles of attack. 
Investigation on the effect of vortex breakdown on the buffeting of 
nearby lifting surfaces, such as tails, was scarce. However, it is 
known that the vortex from the strake (or leading-edge extension, 
LEX) may reduce the buffet intensity on the wing before it bursts 
(p. 109, Ref. 7). 

In the present study, the main objective is to predict 
buffeting on vertical tails induced by LEX vortex bursting. 
Fundamental equations for structural response will first be 
derived. Existing theoretical methods for buffet prediction will be 
reviewed. The present method and some numerical results will then 
be presented. 
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2 . THEORETICAL DEVELOPMENT 
2.1 Formulation of Equations 


Structural Equations of Motion : 

Let the structural displacement, Z a (x, y, t) , be expressed in 

terms of normal mode shapes, <j> n (x, y) . Then 

N 

Z (x, y, t) = l q (t)<j> (x, y) (1) 

a n n 

n=1 

where q n (t) is the so-called generalized coordinates. It can be 
shown that the structural equations of motion in forced oscillation 
in generalized coordinates can be written as (Ref. 13, pp. 131-139, 
or Ref. 14, Chapter 10): 

"nS, + M n“n\ = 2„ = //lP E + P„H„(5. nldCdn (2) 


where 

= //<{>^mdxdy, the generalized mass 

m(x, y) = mass per unit area 

aj n = frequency of the n^ normal mode 

Q n = the generalized force. 

The generalized force consists of two terms, one being the 
externally applied force (i.e., the p_-term) and the other being the 

Jj 

force due to structural motion (i.e., the p M «term) . 

The p M -term can be further decomposed in terms of the 
generalized coordinates as 


« q j 

> M = l APjU, y; 0 ), M J — 

j=1 o 


(3) 
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where Ap^ is the lifting pressure at point <x, y) on the wing caused 
by the motion of the j ** 1 normal mode and b Q is the reference length, 
e.g. the root semichord. it follows that 


N q 


n q, 


Q Mn = ^ b 1 // A Pj<t> n d 5 dn = I ^ If AC ^d 
j =1 o S w j =1 o 


= bjv 


N 


o~ b . ^ A nj^j 
o 3=1 


(4) 


where A^ is the generalized aerodynamic force matrix and is defined 


as 


A . = // AC d> dldn 
ni J J P. n 
: 


(5) 


V 


ti = An 


( 6 ) 


In Equation (4), is the dynamic pressure (= pV /2) . 


Equation (2) can now be written as 
2 N 


M n q + M n°n q n " Aq » ^ = ^ P (5,n,t)* n (5, rOd^dn 

j = 1 


= Q E (t) 
n 


(7) 


In the above derivation, neither structural nor viscous dampings 

have been included. To include the former, is usually replaced 

2 

with oj^ ( 1 + ig^) ' w ^ ere 9 n is the structural damping coefficient 

for the n^ mode and is usually taken to be 0.03 if not known 

experimentally. To account for the latter, 2^M a) q is added to 

n n n 

the equation with £ being the damping ratio. Equation (7) becomes 
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( 8 ) 






N 


M Q + 2?M 0 ) qf + M U) 2 (1 + ig )q n - I A q. = p E (t) 

n ti n n ti nn nn 00 ■ „ ni 1 n 

3*1 


Structural Response to Random Excitation ; 

If the excitation force p n E (t) is random, it may be represented 
in a Fourier integral (Chapter 14, Ref. 15), 


Q E (t) = / Q E ( iw)e iu>t du) 
n 4 n 


(9) 


where 


Q n E (iio) = lim ^ / Q n E (t)e _la)t dt 


( 10 ) 




The displacement q n (t) will also vary randomly, so that a Fourier 
integral representation is appropriate. 


oo , 

q (t) = / q ( ico)e 10Jt da> 
n 4 n 


(id 


Substituting Equations (9) and (11) into Equation (8), and requiring 
the relation to be valid for all t, it is obtained that 


— 2 2 r “ 

q (iu>)[-M a) + 2i£u> u) + M a) ( 1 + ig ) ] - £q ) A .q. 
^n n n nn n ® n] ] 

1 = 1 

= p n E (iu)) 


or, 


N 


Y l[-M u + 2iM a) + M ii) ( 1 +ig ) 1 <5 . - £q A . }q . 

^ 1 n n n n n n nj nj 


= Q E (io)) 
n 


n — 1 , . . • ,N 


( 12 ) 
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Let 


• vT *' V'f •' 


2 2 

2 . (u)) = [-M to + 2iM £ 0 ) 0 ) + M 0 ) (1 + ig )]{ , - tq A , 

nj n n n n n n n] » nj 

(13) 

Note that Z n j ( 0 )) is called the complex impedance of the system? and 
its inverse, Z*” 1 , is the so-called structural transfer function. 

To describe quantitatively a random response in a meaningful 
manner, statistical methods must be used. The most important 
quantity for this purpose is the mean square value. It is defined 
for a random function F(t) as (Ref. 15) 


■ i A w . 

F 2 (t) = lim -L J F 2 (t)dt = lim -L / F(t) / f ( iu))e 1U)t dtodt 

X-^oo —X T+°° — T — oo 


= lim / f(iai) / F( t)e itot dtdaj 

Orp •' * 


X+° 


2T 


-T 


00 oo 2 

= lim L- / f(ia))2irC (iu)do) = / lim — duj (14) 

2T J J T 


X-Voo 


— OO X“^°° 


where f is the complex conjugate of f. Define 


S( to) = lim - 
X-K» 


irf(i(o)f (ito) 


(15) 


Equation (14) becomes 


oo 

F 2 ( t) = / S(to)dto (16) 

— OO 

in case the random function depends also on space coordinates, 
the definition of S(io) must be modified. For the generalized force 
of the n th mode, Q n E (t), it is defined as (Equation 7) 
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Q E (t) = / p (r, t)(J> ( r )dA 
n E n 


( 17 ) 


where space coordinates (£, ri) are now represented by r. The 
Fourier spectrum of Q n E is 


^ B (i») = // P E (iu)) <J> n (r)dA 


( 18 ) 


The power spectrum of the n th generalized force is given by 


S (w) = lim 5 E (iw)Q E (ito) 

n T ti n 

T+°° 


- li ® § f ^ E ( i “ ) * n (r 1 ,dA 1 ^ P *< i “>* n < r 2 ,dA 2 

T-v» A A * 


T+oo AA 


' I // *„ (r , 'W 'h F C P E <r i ' *1 ,P E ( V V 

4tt -T -T 


ico( t -t ) 

• e dt dt dA. dA„ 

12 12 


( 19 ) 


Let t 2 - t 1 = t. Equation (19) can be written as 


S (u>) = // 4> (r ) 4> (r )lim - — / / p (r , t )p (r , t + t) 
n ii n 1 T n 2 _ T 4i „ „ E 1 1 E 2 

AA -T -T 


T T 


- 10)1 

• e dt dxdA dA 
1 1 2 


where 


■ [[ WW 1 iiW V T)e ' 1 “ TdTda , da : 

AA 


( 20 ) 


, ,2 <r r r i T) ' u "-k / W WVS * ,)d S 12,1 


T-^-OO -T 

is known as the space-time correlation or cross correlation 
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function. Define 


S 12 <V V “)-■??/ r , 2<V V T|e ' 1 “ TdT (22 > 

.00 

the cross power spectral density of pressures at r 1 and r 2 « It 
follows that 

V“> - I W / *n (r 2 ,S 12 (r r r 2 “ )d V A l 1231 

A A 

S(w) or S n ((o) is known as the power spectral density. This is 
because if F(t) were a current, the power developed by this current 
as it passed through a resistance of one ohm would be F 2 (t). 

Returning to calculation of the total response, Equations (11) 

- (13) show that the amplitude of the motion in the n** 1 mode is 

00 . 

q n (t) = / ([ZU)]' 1 {Q E (iu>) }) n e lt0t du) (24) 

— OO 

The total displacement is therefore 

N 00 « __ • t 

Z (x, y, t) = y / r [ Z ( to) ] ' {q ( ioo) } ) <p (x, y)e doo (25) 

a n n 

n=l 

from which the Fourier spectrum of the total displacement can be 
identified as 


l ([Z( u) ] 1 {Q E (iio) } ) n 4> n (x, y) 
n=1 


(26) 


and the corresponding power spectrum is 


s „‘“> 


N N 

lim | { l ( [ Z( co) ] ~ {q (iu>)}) n <j> n (x, y) }{ £ ([Z(u)l 
T^°° n=1 n=1 


{q E (±o>) } ) n * <t> n (x, y)} 


(27) 
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Once S w ( uj) is known, the mean square value of displacement can 


be obtained as 


Z 2 = J S ( a>)du) 
a J w 


(28) 


Responses in accelerations, loads, moments, and stresses, etc,, 
can be similarly formulated. 

Equation (27) is difficult to simplify because of mode coupling 
through the generalized aerodynamic force matrix, A^j . If the 
aerodynamic force due to structural motion is ignored, or A^ ® 0 if 
n ^ j, then Equation (27) can be further simplified. Let 


Z ( (d) = M [- 0 ) 2 + 2i£0) uj + <o 2 ( 1 + ig ) ] - £q A 
nn n n n n °° nn 


(29) 


Equation (27) can be rewritten as 

N p E (ito)<b (x, y) N 5 E (iu))<f> (x, y) 

S W U) - { l -2 }{ l ^ 

T>°° n=1 Z ( w) n=1 Z ( u>) 


} 


nn 


nn 


(30) 


After multiplying this out, it can be obtained that 

_ N Q _ E Q®* 2 ( x , y) • 

S W < W> = lim T 1 i * + 

T>°° n=l Z Z 

nn nn 

- E- *E 
N N Q. Q n 

+ l l ~ y t ] 

j=1 £=1 Z . . Z„ . J 

Ml ££ 

N <f> 2 (x, y) 

' \ 17 T 7 T l VV l VVW r 2 ' “ )dA 2 dA 1 * 

n=1 Z ( to) A A 

nn 


N N ^(x.yJUx, y) 

*11 ~ / ♦ j < r 1 ) / ^ (r 2 )S l 2 (r r V u)dA 2 dA 1 


j=1 £=1 Z . . (tb)Z (to) A 


(31) 


9 


The first series of Equation (31) represents the sum of the spectra 
of the responses in individual modes. The second series represents 
the correlation between the responses in different modes. The 
second series can be ignored if only two or three modes are present 
and their natural frequencies are widely separated (Ref. 15). 

in Reference 16, the cross power spectral density was specified 
in exponential functions with coefficients determined by experiment. 

2.2 Existing Theoretica 1 Methods for Buffet Predictio n 

All existing theoretical methods require some types of 
experimental data to work with. Sophistication of these required 
data distinguishes one method from the other. 

Cunningham and Benepe (Refs. 17 and 18): 

Pressure power spectral densities are first converted into 
pressure distributions over the wing for each frequency. The 
doublet lattice method ( DLM) is then used to calculate induced 
pressures on the tail due to downwash produced by the wing buffet 
pressures. The wing and tail pressures are used in the DLM to 
calculate the generalized aerodynamic forces. The whole equation 
(12) is used without simplification. The calculation is similar to 
that for gust response. 

B. H . K. Lee (Ref. 19): 

Again, Equation (12) is used. However, the cross correlation 
function S 12 in Equation (23) is either taken to be constant over an 
aerodynamic panel or asumed to vary exponentially in space. 
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Mullans and Lemley (Ref. 20): 





The fluctuating pressure on a rigid model is again used to 
calculate the generalized aerodynamic forces. However , the 
aerodynamic forces due to wing vibration, i.e. A^-terms, cure 
ignored. 

J. G* Jones (Refs. 21 and 22): 

It is assumed that each mode behaves as a single-degree-of- 
freedom system: ^ 

"A + + vA - <A> (32> 

The aerodynamic forces due to wing motion are ignored. Applying the 
Fourier transform to Equation (32), it is obtained that 


M (-w 2 + 2i?u> uj + aj 2 )q =p E (iw) 
n n n A A 


(33) 


Using the definition of power spectral density, Equation (15), the 
power spectral density of the response can be obtained: 


S (to) = — T ^ 

q .. 2. 2 


s p (a)) 


^ 2 2 2 
M *(-u>~ + 2i£w a) + to ) ( — co - 2i£u) 0 ) + a) ) 
n n n n n 


(34) 


The mean square value of q n is therefore 


q ( t ) = / S ( w)dw 
n ' q 


- / 


S- (») 
n 


2 * 

M H ( 0)) H (U>) 
n n n 


doj 


(35) 


where 


2 2 
H (w) = -u> + 2i?U) a) + 0) 

n n n 


(36) 
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The main contribution to the value of the integral in Equation (35) 

comes from the peak response at id = u , if s-(to) is assumed not to 

n Q 

vary appreciably in the neighborhood of u) n , it can be factored out 
of the integral in Equation (35) and the result integrated 
analytically based on the residue theorem in the theory of a complex 
variable. Results are available in Reference 23 (p. 218). 

Therefore, Equation (35) can be reduced to 


— s e 'V 

(t > ‘iA—i 


(37) 


2 ”2 

Instead of q^ (t), Jones determined q^ , the mean square 


acceleration. Note that the Fourier transform of q is 

n 


q = (ico) 2 q 
n n 


(38) 


Therefore , 


s q .. 2 


(ito)*S- (to) 

% 


M H (w)H (a)) 
n n n 


(39) 


•• 2 

The result for q is (Ref. 21) 
n 


• * * co 

2 1 n _ . . 

q r> “ P T~ S n < *0 

n 8 M 2 C Q n n 

n 


(40) 


Let 


S- (co) ( q<o S) 2 < 

v n 

where is the freestream dynamic pressure and E 2 is a 
nondimensional aerodynamic excitation parameter. It follows that 
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) 


( 42 ) 


« {' /2 q 2 

E - 2/2 <JL>' /2 < »„ =- 

- 0 ) S q 

c n ® 


i c “> , ,, q s 
J_ ( _il) 1 /2 e ~ 


" 2/2 VC 

To use Equation (43), the damping ratio (c;) is needed* It 


(43) 


consists of both the structural damping ( f g ) and the aerodynamic 
damping (£_). The latter arises from the effective angle of attack 

d 

due to wing vibration and is given by 


2M ^ 0 ) q = 2q SK — 
n a nm 00 V 


(44) 


where K, the aerodynamic damping parameter, is a nondimensional 
parameter depending on the mode shape, the wing planform, and the 
sectional lift-curve slope* Equation (44) is assumed applicable to 
both attached and separated flows* It follows from Equation (44) 
that 


” M a) V 
n n 


K 


M a) VC 
n n ^a 

qTs”” 


(45) 


(46) 


In Jones 1 method, both E and K are assumed to be independent of 
the scale effect. in other words, their values determined from 
model test can be applied to full-scale airplanes. practical 
procedures of applying this method were discussed by Butler and 
Spavins in Reference 24. They are as follows: 

(a) Determine modal frequency u) n , the mode shape, generalized 
mass v and structural damping from wind-off resonance 


tests on model and aircraft. Note that the relevant model 


mode shape must be approximately correct. 

(b) Measure rms acceleration or bending moment c~~ at a point 

B 

on the wing, the total damping £, flow velocity v r and 
dynamic pressure q^, at a given Mach number and angle of 
attack in wind-tunnel tests. 

T\ 

(c) Relate to q^ in generalized coordinates using the mode 
shape (see Section 2.3). 

(d) Calculate E from Equation (42). 

(e) Calculate k from Equation (46). 

(f) Calculate total damping of aircraft by adding calculated 
£ a from Equation (45) to the measured 

(g) predict rms acceleration or bending moment at a point on 
the aircraft wing from Equation (43) using the measured 
aircraft mode shape. 

Mabey's Method (Refs. 25 and 26): 

This method was developed to determine qualitatively the flight 
conditions for light, moderate and heavy buffeting for the full- 
scale aircraft from measurement of wing root bending moment of a 
conventional wind-tunnel model. It is assumed that the wing 
responds to buffeting pressures in somewhat the same way as to the 
wind-tunnel turbulence at the wing fundamental frequency. 

Let the tunnel unsteadiness /nF(n) be defined so that the total 
rms pressure fluctuation coefficient is given by 
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2 


(47) 


P 

q 


00 


oo 



(nF( n) ] 



where 

n = f w/V 
o 

w = tunnel width 

f = wing fundamental bending frequency in cycles 
per second 

V = frees tream velocity. 

Define 


C BB^ M ' = wing-root strain signal/q^ 


(48) 


Before the onset of flow separation on the model , 


C 


BB 


(M, 


a) has been 


shown experimentally to be constant equal toC (M, a = 0). This is 

BB 

the portion of the model response caused by the tunnel unsteadiness /nF(n) 


. Assume that 


C bb (M, a = 0°) = Kg/nFfn) ( 49 ) 

where K g is a scaling factor. Then 

C ’(M, a = 0°) = — C (M, a = 0°) = /nF( n) 

&& K BB 

B 

(50) 

Beyond buffet onset, Cg B (M,a) is increased due to wing buffet 
pressures. Let 

C bb "(M, a) = [C bb '(M, a) 2 - C bb '(M, a = 0°) 2 ] 1/2 (51) 

The angle of attack at which C BB " (M, a) first differs from zero is 
the buffet-onset angle. From correlations on nine models of fighter 
aircraft, the following buffeting criteria were suggested: 
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Buffet onset 


0 


Light buffeting 
Moderate buffeting 
Heavy buffeting 

Note that in using this method, the total damping of the wing 
fundamental mode should be relatively constant, independent of wind 
velocity and density. This is true if models with solid wings of 
steel or light alloy are used, because in this case the structural 
damping will predominate. No mass, stiffness (or a) o ) and damping 
for both models and aicraft are needed. It is useful during 
comparative tests for projects with alternative wing designs. 

Thomas* Method (Ref. 27): 

At transonic speeds, buffeting is closely connected with flow 
separation due to shock-boundary layer interaction and shock 
oscillations. Using conventional boundary layer methods, the 
development of boundary layer on airfoils at transonic speeds can be 
calculated. By comparing calculations with experimental results, 
Thomas postulated that buffet onset started if the point of rear 
separation coming from the trailing edge reached 90% of the airfoil 
chord. 

Redeker (Ref. 28) extended this method to infinite yawed wings 
by using the pressure distribution on a section normal to the 
leading edge and applying a three-dimensional compressible boundary 
layer method. 


c BB " = °* 004 
c bb" = °* 008 

C BB n = °* 016 
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Further extension of Thomas* method to finite wings was made by 
Proksch (Ref. 28). A buffeting coefficient (C Bi ) is defined which 
is directly related to the rms value of the wing root bending 
moment. It is assumed that the fluctuations of the wing root 
bending moment are proportional to the integral evaluated along the 
wing span of the product of local lift fluctuations and the distance 
from the wing root (n - n^) • A further assumption is that the local 
lift oscillations caused by flow separation are proportional to 
length £ g (n) of the separated flow at a spanwise station of the 
wing, it follows that 

1 —5- 

c Bi = / — <n - ri R )dn ~ /c B (52) 

. n R C 

2.3 The Present Proposed Method 

Theoretically, it is possible to use Equations (23) - (28) to 
calculate buffet response in the most general way. However, it 
would be an expensive undertaking because extensive fluctuating 
pressure measurement on empennage must be made. In addition, these 
fluctuating pressures are configuration dependent and vary with 
flight conditions. Therefore, a method similar to Jones* in concept 
is proposed. 

In developing the proposed method, the following steps are 
needed. 

(a) Buffeting vortex strength in the burst region must be 
known. It is known that steady vortex strength from a 
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slender wing or LEX can be estimated by the method of 

suction analogy (Ref. 29). Similarly, buffeting vortex 

strength can also be estimated if buffeting normal force 

data on slender wings are available. This is because any 

buffeting on slender wings can be assumed to be caused by 

the leading edge vortex. A limited amount of such data 

was published in References 4 and 9. 

Let c e be the sectional suction coefficient. Based 
s 

on the suction analogy, the vortex lift is proportional to 
c g . The vortex lift can also be expressed in terms of the 
vortex strength r through Kutta-joukowski theorem as 

pV 2 c cdy = pIV d£ (53) 

2 00 s £e 

where w^ e is the normal velocity at the leading edge and 
d£ is the vortex length along the leading edge. It 


follows that 


c c 
s 


f d£ - i tht dy 

00 £e 00 


and 


r t -/f 


b/2 c c 
s 


w„ A 

£e o 


dy 


The average strength per unit length is 


1 b/2 c c 

f t ‘ V S «e ’ 2S - I » 7-^- ^ 

£e o 


■* e /v - 


(54) 


where S^ e is the length of the leading edge. The unsteady 
aerodynamics program of Reference 30 was revised to 
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calculate f^. In the calculation, the buffeting normal 
force is obtained by assuming a vertical oscillation of 
constant amplitude over the region of predicted vortex 
breakdown. The latter was calculated by a semi-empirical . 
method to interpolate or extrapolate experimental data 
(Ref. 31). The amplitude was adjusted to match the 
experimental data on fluctuating normal force coefficients 
in Reference 9. Unfortunately, only data at a low 
frequency for some delta wings were measured in Reference 
9. On the other hand, the power spectrum over a range of 
frequencies at the vortex-breakdown angle of attack for 
the BAC 221 configuration is available (Fig. 24 of Ref. 

4). Unless additional data are available in the future, 
for the present purpose the low-frequency data of 
Reference 9 will be used to derive the buffeting strength 
for a range of angles of attack. At other frequencies, 
the strength will be multiplied by a ratio obtained from 
data for the BAC 221 in Reference 4. 

(b) The root bending moment can be calculated as 

b/2 N .. 

M Q ( t) = / U E (y, t) + * M (y» t) - ( l q n (t)<f> n (y) )m(y)]ydy 
o n=1 

( 55 ) 

where is the sectional lift due to external forces, 
is the sectional lift due to structural motion and the 
last term is the inertial forces. For a rigid wing, the 
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last two terms can be ignored. In Jones* analysis, was 
also ignored. Using the notation of Equation (4), can 
be written as 

N N q. 

*M (yr = q oc I I j^- / AC (x, y)(4> n (x, y)mdx (56) 

n=l j=1 o p j 


Let the Fourier transform of £ £ (y, t) be written as 


N 


* E (y, t) = l (y* i( 0 )Q n E (iw) 
n=1 n 


(57) 


where & is the sectional lift due to a unit generalized 
E 

n 

force in the n** 1 mode. Applying the Fourier transform to 
Equation (55), it is obtained that 


b/2 N £ N 2 

M o (i 0)) = / [ ll (y, iu) + — + l Tr |uT m<fr n (y))Q„ ydy 

o n=1 n O n=l n 

*n 


n 


N 


= y [BM (iw) + BM (iu) + — 
E M M 

n=1 


0) 


n 


n 


H (u>) 
n n 


“ E 

/m<}> n (y)y dy]Q n 




(58) 


where Equation (33) has been used. H BM (co) is the bending 
moment transfer function and is defined as 


N 


H (u>) 
BM 


= y ( BM (ito) + BM (ico) + — 
E M M 


0) 


b/2 


n=1 


H <«) / "*n ly)y dyl 
n n o 


( 59 ) 


The power spectral density of M Q (t) is therefore 




where S- is the power spectral density of the buffeting 
excitation. For a rigid wing. Equation (60) can be 
simplified to 


(w) (61) 

in applications, BM^, will be calculated by assuming a 
unit buffeting excitation over the region of vortex 
breakdown at a range of frequencies. The mean square 
value of bending moment is then given by 


V - \ 


M = / S (u>)du) = 2 / S_(io)dio (62) 

o 1 BM J BM 

— oo O 

which is to be integrated numerically. 

(c) Since only total force power spectrum, instead of pressure 
power spectrum, will be used, it is assumed that the 
pressure fluctuations at every point on the wing are 
perfectly correlated in space and are in phase. Based on 
this assumption, Mabey and Butler showed that the total 
force power spectral density was proportional to the 
pressure power spectral density (Ref. 11). The results 
from this were shown to be reasonably accurate. 

In the present application to empennge buffeting due 
to a LEX vortex, those unsteady buffeting vortices, once 
generated, will be convected downstream in accordance with 
the general principle of unsteady aerodynamics. 

(d) With the power spectral density of buffeting vortex 
strength determined at a given flight condition, 
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fluctuating normal velocity will be induced on the 
empennage. By satisfying the usual flow tangency 
condition, buffeting pressure spectral density on the 
empennage can be calculated. From the buffeting pressure 
spectral density, the power spectrum of bending moment or 
other aerodynamic characteristics can be determined. The 
root mean square values of root bending moment are 
calculated by using Equation (62). 

(e) Similar to Jones 1 method, the calculation of buffet 
response requires structural data, such as mode shapes, 
generalized mass, and damping ratio, 

(f) in applications to empennage buffeting, the locations of 
LEX vortex bursting will be based on experimental data. 


22 


3 • SOME NUMERICAL RESULTS AND DISCUSSIONS 

In Reference 9, data for a cambered and a flat model of 70- 
degree delta wing are available at a frequency parameter n of 0*05, 
where 

n = fc/v (63) 

oo 

and f is the frequency in cycles per second. This is converted into 
the conventional reduced frequency k by multiplying by 2 tt. Before 
vortex breakdown, the normal force fluctuation is assumed to be 
caused by tunnel flow unsteadiness. The resulting amplitudes of 
vertical oscillation (or buffeting excitation) in the vortex- 
breakdown region needed to produce the experimental mean square 
values of normal force coefficients are shown in Figure 1 . At each 
angle of attack, a buffeting vortex strength f can be calculated 
from Equation (54). The same expression is used to calculate the 
vortex strength F in steady flow using the steady-flow c g . Now, if 
Figure 1 is replotted in terms of the ratio of buffeting to steady 
vortex strengths, R b/ S : 

Vs ■ V T s <55 > 

those two curves in Figure 1 tend to collapse into one as shown in 
Figure 2. From Figure 2, it can be concluded that the buffeting 
vortex strength is a function of steady-flow vortex strength and Aa, 
where Aa is the incremental angle beyond that of vortex breakdown 
which occurs at the trailing edge. 


23 


The meth6d is now used to analyze a cambered delta wing of 65- 
degree sweep as shown in Figure 3. To find the buffeting 
characteristics for a flat 65-degree delta wing, it is assumed that 
the buffeting excitation (AMPLG) for the latter is equal to that for 
a cambered 65-degree delta wing if R^g is the same. Therefore , at 
a given Aa, is obtained from Figure 3. Using this R b / s , AMPLG 

can be determined from Figure 4* Note that Figure 4 was constructed 
from the experimental data for a cambered wing. The resulting 
buffeting excitation amplitudes for a flat delta wing are plotted in 
Figure 5. This type of data extrapolation must be used with 
caution. This is done here mainly because test data for fluctuating 
normal force coefficients of the flat 65-degree delta wing are not 
available. However, test data for the rms values of root bending 
moment are available for comparison with predicted results. 

The buffeting force characteristics for both 70-degree and 65- 
degree delta wings are compared in Figure 6. The difference between 
these two curves arises mainly from differences in forward 
progression rate of breakdown points. 

To check the theory, test data of Reference 32 for the root 
bending moment will be used. Static bending moment coefficients 
based on c are presented in Figure 7. Calculated results from 
Reference 29 are also presented for comparison. It is seen that at 
high angles of attack, the theory overpredicts the root bending 
moment. This is because the theory does not account for the inboard 
movement of vortex flow as the angle of attack is increased. An 
investigation to correct this discrepancy is under way. 
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To calculate the dynamic response of a rigid wing, Equations 


(61) and (62) will be used. To calculate the transfer function 
(BM e ) for the root bending moment, a unit amplitude of vertical 
excitation is prescribed over the region of vortex breakdown at each 
frequency. Some results are presented in Figure 8. The 
corresponding power spectral densities for the excitation are 
obtained by multiplying the values in Figure 5 (for a low frequency 
only) by a ratio obtained from Figure 24 of Reference 4 for other 
frequencies. The results are shown in Figure 9. Equation (62) is 
then integrated by the trapezoidal rule to produce the mean square 
values of root bending moment. The rms values are presented in 
Figure 10* Note that experimental data were obtained at resonant 
frequencies of the fundamental bending mode. Since the spectral 
density is higher at higher frequencies (Fig. 9), the calculated 
response of a rigid wing tends to be similar to the test data at a 
high frequency, although the magnitudes are underpredicted. It is 
expected that the prediction can be improved if the structural 
flexibility is accounted for. 
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4 « WORK UNDER WAY 


All geometric data for the F-18 have now been acquired and 
computerized in the input format for the VORSTAB code (Ref. 29). 

Wing sectional aerodynamic characteristics are being calculated with 
the Eppler's code (Ref. 33) to form a part of the input. At a given 
angle of attack, the wake and the LEX vortices will be allowed to 
deform until an equilibrium position is reached. The resulting 
position of the LEX vortices relative to the vertical fins of F-18 
will then be input to the unsteady aerodynamics program of Reference 
30 to calculate the fin buffeting. For this purpose, the program of 
Reference 30 must be modified to include the effect of large 
geometric dihedral angles. Again, only the response of rigid fins 
will be calculated in the current plan. 
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Appendix B: 

Water Tunnel Testing of an F-18 Model 
by 

William H. Wentz, Jr., 

Wichita State University 


DISK: S5 FILE : FA18/JA6 

SUBJECT: PROGRESS REPORT, F/A-18 FIN BUFFET PROJECT 
To: Eddie Lan, KU Date: 20 Jan 1986 

From: Bill Wentz, WSU 


STATUS OF WORK STATEMENT ITEMS: 

The work statement items for the WSU portion of this 
project are quoted from the proposal, followed by the 
status for each item. 

1. "Participate in detailed water tunnel test planning 
with NASA-Dryden and Navy personnel."- Planning and tests 
were completed during December 1985. 

2. "Assist in the interpretation of available 
literature on F-18 and similar aircraft pertinent to the 
fin buffet problem."- This activity is to be in 
coordination with Professor Lan. No activity to date on 
this task. 

3. "Review data obtained from water tunnel tests. "- 
This has been a continuing process since completion of the 
tests. Key results are included in later sections of this 
report. At the present time, no additional water tunnel 
tests are recommended. Possible recommendations for 
addtitional water tunnel tests, and recommendations for the 
flight tests will follow completion of analysis of water 
tunnel test data. 

4. "Contribute to the preparation of a technical 
report summarizing the results of this project."- Water 
tunnel test results are being prepared in form appropriate 
for the final report. 


FACILITY AND MODELS: 

Water Tunnel tests were conducted in the NASA-Dryden 
flow facility during the period 13-17 Dec 1985. Model 
configurations were as follows: 

Basic F/A-18. - This model is essentially the complete 
aircraft configuration with missiles removed. Wing leading 
edge flaps were deflected 34° and trailing edge flaps were 
un-def lected. This configuration is consistent with flight 
operations at angles of attack above 25°. 

F/A-18 without wings. - This model was used to 
evaluate the interference effects of the wing and leading 
edge extension (LEX) flow fields. 
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F/A-18 without fins. - This model was used to evaluate 
the possibility that the fin "blockage" might generate an 
adverse pressure field of sufficient strength to cause 
premature bursting of the leading edge vortices. 

F/A-18 without LEX's. - The purpose of this model was 
to identify the role and interaction of forebody vortices, 
and to ascertain possible wing or forebody vortex 
interactions with the fin. 

All models were constructed from 1/48 scale 
commercially available Monogram kits, fitted with 
hypodermic tubing for introduction of flow tracers. The 
models were also equipped with engine exhaust tubing 
connected internally to the engine inlets so that engine 
inlet flow was simulated. 


TEST CONDITIONS: 

For all but a few tests, inlet flow was established to 
provide an inlet capture area ratio of one. Most tests were 
conducted at a tunnel speed of 0.25 ft/sec, which 
corresponds to a Reynolds number of 4,000 based on (1/48 
scale) model wing mean aerodynamic chord of 0.24 ft. Mach 
number for this speed is 5 x 10~ . 

Angle of attack was varied from 0° to 40°, in 5° 
increments. At 40°, the model nose was nearly in contact 
with the test section wall, so higher angles could not be 
tested without use of an offset sting. Video and still 
photos were obtained from top and side views in separate 
runs. One series of runs were made with the basic 
configuration with 5° sideslip. 


INSTRUMENTATION: 

Instrumentation consisted of video cassette recording 
equipment, and camera for still photos. In addition, the 
fin of the basic model was fitted with a strain-gage and 2 
different types of surface hot-film anemometers (Disa and 
Micro-Measurements). These instruments were intended to 
detect unsteady flow over the fin, for correlation of fin 
buffet with flight test data. An oscilloscope was available 
for monitoring of strain-gage or hot-film output signals, 
and a modal analyzer was utilized to perform frequency 
analysis of the dynamic signals. 

The Disa hot-film anemometer and the strain-gage 
provided only very low-level signals, and did not provide 
consistent-results which could be distinguished from random 
noise. The MM hot film gage, however, produced a signal 
which displayed characteristics which changed in a consistent 
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manner with angle of attack. Therefore only the data from 
this gage was utilized for spectral analysis. 


RESULTS OF FLOW VISUALIZATION: 

BASIC MODEL - The flow video and still pictures show a 
consistent and repeatable pattern for the vortex flow of 
this aircraft configuration. As angle of attack is 
increased from 0°, increasingly stronger vortices form 
along the LEX's. These vortices flow aft above the 
horizontal tail surfaces but beneath the fins for angle 
below 20°. At 20°, busting occurs aft of the wing trailing 
edge and outboard and beneath the fin. At 25° angle of 
attack, the LEX vortex burst point is located inboard, with 
the primary axis of rotation nearly coincident with the fin 
leading edge. The burst point is slightly forward from the 
fin leading-edge at this angle of attack. As angle of 
attack is increased further, the burst point progressively 
moves forward. Vortex burst locations are shown in figures 
1 and 2. 

MODEL WITHOUT FINS - The absence of the fins had 
little effect on the flow field. Vortex locations and burst 
position were essentially unchanged from the basic model. 

MODEL WITHOUT WINGS - In this configuration, LEX 
vortices formed in much the same manner as for the full 
model. As angle of attack was increased, however, the 
vortices were located at more inboard location than the 
basic model, and they remained intact, without bursting, up 
to 30° angle of attack. This test series clearly indicates 
that the wing pressure field has a dominant role in the 
vortex bursting process. The adverse pressure gradient 
field associated with the portion of the wing aft of the 
leading edge flap hingeline is evidently a dominant factor 
producing vortex bursting. 

MODEL WITHOUT LEX'S - With this model, no clear 
picture of vortices impinging on the fins was observed. 
Comparison of this configuration with the basic model 
reveals the strong role of the LEX's in producing the 
vortices which impinge on the fins. 


RESULTS OF HOT-FILM SIGNALS: (BASIC MODEL ONLY) 

The modal analyzer was utilized to obtain Psd data for 
each angle of attack from 0° to 4 0°. Results of these 
studies are summarized as follows: 

(1) For angles of attack of 0° to 20°, no 
dominant frequency was observed. 

(2) For angles of attack of 25°, 30°, 35° and 
40°, dominant frequencies were discernable. 
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For 30° and 35°, runs were also made with tunnel 
speeds of approximately two- and three-times the nominal 
value. Results were plotted as frequency versus velocity in 
figures 3 through 6. These results show that frequency 
increases linearly with tunnel speed, resulting in a 
constant Strouhal number. Further, the Strouhal number 
associated with the vortex bursting is essentially 
independent of angle of attack. The observed Strouhal 
number is approximately 0.7 for all cases. 

All test conditions scheduled have been run, and all video 
cassettes, photos and modal analyzer plots have been 
provided to WSU. 

CURRENT ACTIVITIES: The NASA-provided video 3/4-inch tape 
cassettes have been transcribed onto 1/2-inch tape 
cassettes for ease of analysis using the WSU stop-action 
VCR unit. Analysis of these tapes is continuing, and WSU 
will extract primary LEX vortex location (vortex core 
location and location of the bursting point) from these 
images for the sideslip and non-standard configuration 
cases. Additional narrative describing the flow phenomena 
and associated hot film measurements will be developed. 
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